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Abstract 

We apply the Keldysh formalism in order to derive a current formula easy to use for a system 
with many sites, one of which is interacting. The main technical challenge is to deal with the 
lesser Green function. It turns out that, in the case of the left-right symmetry, the knowledge of 
the lesser Green function is not necessary and an exact current formula can be expressed in terms 
of retarded Green functions only. The application is done for a triangular interferometer which 
gives a good account of the Fano-Kondo effect. It is found that the interference effects, in the 
context of Kondo correlations, give rise to a point in the parameters space where the conductance 
is temperature-independent. We include a comparison with the results from the Ng's ansatz, which 
are less accurate, but can be used also in the absence of the above mentioned symmetry. 

PACS numbers: 73.23.-b,73.63.Kv,85.35.Ds 



1 



I. INTRODUCTION 



The transport through mesoscopic systems has been very much discussed in the last years 
because of promising technological applications, but also for revealing interesting conceptual 
aspects. The simultaneous consideration of both interaction and interference effects is a 
nowadays topic in mesoscopic physics and many efforts have been done to overcome the 
specific theoretical difficulties. One of the main tools is the Keldysh transport formalism, 
that requires the knowledge of the retarded and lesser Green functions (e.g. [sl, I^, 0, [s]). 
The general expression of the current through the lead indexed by a is : 

^ = 7^ y du tV^{Uuj)[Gl^{u) - + G<^{uj)} , (1) 

where fa is the Fermi distribution function in the a electrode and Fq, = vrr^p, with the usual 
notations. 

Most of the papers discuss the Single Impurity Anderson Model (SIAM) , in which case a 
simplified version of the formalism can be employed. In this model, the two leads are con- 
nected to the same site, which mimics the quantum dot, and in this situation the knowledge 
of the retarded Green function is sufficient for describing the transport properties (as the 
lesser Green function is eliminated by current symmetrization)[l|. Only a few papers go 
beyond this model by introducing in the Hamiltonian a term which produces a short-cut of 

riQ 

the impurity, mimicking a mesoscopic ring [6|, [Zl]. In the non-interacting case, eq.(l) can be 
shown to give the Landauer-Biittiker result, that was widely employed. 

The implementation of the Keldysh formalism for the general case of a many-site, multiple 
connected system with interaction is not trivial. Second order perturbative calculations 
have been applied for the case of the electron-phonon interaction j^, 4| or in the context 
of dephasing jsl, but such perturbative approaches cannot be applied in the case of strong 
electronic correlations. Recently, Entin-Wohlman et al. [s^ and Kashcheyevs et al. 9] 
developed an equation-of-motion solution for G^' in any complex geometry described by a 
tight-binding model, but not for . The calculation of G^ in the strong interaction regime 
is a more difficult task and requires supplementary decoupling approximations j^. [lo|. 

A solution to this problem was also suggested ll|, [l2| by diagonalizing the interacting 



Hamiltonian in the slave-boson representation. Then the knowledge of the retarded Green 
function is sufficient again, but the approach restricts the range of validity of the transport 
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calculation as the mixed valence regime is not correctly described. Another way to get 
rid of was proposed in [13] by imposing stationarity conditions on the mean values of 
some non-hermitic operators, the meaning of which is not obvious to us (especially under 
non-equilibrium conditions). 

In this paper we propose two different approaches. One is to approximate the lesser Green 
functions by using the Ng's ansatz [14]. Another possibility is to write down the system of 
equations of motion for Gfj (where i,j index any site of the system) and express the current 
(1) in terms of and (the index 'd' means the interacting site). It will be shown that, 
for the particular case of left-right and time-reversal symmetry, can be eliminated by 
symmetrization of the current formula. The result is better in this case than by using Ng's 
ansatz, since no approximations are introduced in the current formula. 

The application is done for the triangle system in the inset of Fig.l which is the simplest 
multiple connected system with two leads. This system is typical for the study of Fano and 
Fano-Kondo effects. 



II. MODEL AND TECHNIQUES. 

Our approach is based on the following Hamiltonian written in a discrete basis : 

k,(T,Q! 

where the mesoscopic system contains an interacting site indexed by d: 

Hmeso = ^ Eial^ttia- + "^{Ujal^aj^ + H.c.) + Hint{a\^, ada) (3) 

ia ija 

and the tunneling term is of the following form 

Ht=J2 (^^ia,a(^c.,a + H.C.), (4) 
k,cr,Cl 

where a is the lead index (in c|^^ ^), but also stands for the site where the lead is attached (in 
a-a,cr)] Ta is the coupling constant. Hint may describe any interaction which implies the site 
"d", such as Hubbard or electron-phonon (photon) interaction, the electrostatic coupling 
to external^detectors, etc. The equations of motion for the lesser Green functions can be 
written as 



J2 (^ik - g:Uk)G< = Tlg<glGt^ + glUdG^^ , i,]^ d. (5) 



In the above equation we use the following notations: = {u — Ei + iTi)~^ and = 2i7ipfa 
is the lesser Green function of the lead a which is coupled to the site "z"; p is the flat band 
density of states. The solution of eq.(5) reads 

G<. = {A-'),, {rlg<glGl^ + ^7l4.G< ) , (6) 

where the summation over "fc" is assumed, and the notation Aik = Sik — g^tik have been 
used. The function G^- in eq.(6) can be obtained in a similar way: 

G< = {{AT'),k{^l9^9tG:,, + gtUkG<,) (7) 

where g'^ = {gD*. 

One notices that Gfj can be expressed in terms of the retarded (advanced) Green func- 
tions and G^^. The calculation of G^^ still remains a problem and several approximate 
solutions can be used, as for instance the Ng's ansatz described later. However, an impor- 
tant simplification occurs in the symmetric two-lead case when, after the symmetrization 
J = {Ja — Jp)/'^, the current becomes independent of : 

J='^ j dujTaifa- ff3)l -lraGl^+ 

(8) 

The above equation represents the main formal result of this paper. The current J was 
expressed solely in terms of different retarded Green functions, for the symmetric case. The 
symmetry is necessary in order to get the same coefficient of both in the expression of 
and Ja, so that is eliminated by the symmetrization (J^ — -^/3)/2 . One can notice that 
the SIAM formula is recovered by the first term in eq.(8). Eq.(8) contains many retarded 
Green functions, but, in fact, all of them can be expressed in terms of G^^, as for instance: 
Gd-y = {^'^)jp9ltdpG''^^. 

Ng's ansatz . The lesser Green function G^ = G^'H^G"' is approximated by assuming 
that = S°^M |l4| where the matrix M is deduced from the relation — = S*" — E*^ 
. The result is : 

S< = S°<(S°"-S°")~i(S"-S"). (9) 
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The intention is again to keep in the final formula only retarded quantities and non- 
interacting functions (trivial to calculate). The non- interacting selfenergies , for the two 
lead system read : 



fa \ 




U ) 






\ 




I r J 



(10) 



It is now straightforward to express the quantities required by the current formula eq.(l) as 

Finally, the current formula becomes : 

Ja = -j-l dcoif^ - ^)r„G;^[(S^' - S")G1^„. (12) 

The same scheme was used in jla . IitI , but for a different problem, namely the spin transport 
through a single-site dot coupled to magnetic leads. 



III. APPLICATION AND DISCUSSIONS. 



The exact formula eq.(8) will be applied to the particular case of a triangular interferom- 
eter with Hubbard interaction Hint = f^'^d^'^dT'^di'^'ii • -'■^ previous section, we reached 
our goal of expressing the current through the retarded Green functions only. One impor- 
tant advantage is that there are already recipes for computing these functions in different 
approximations. We shall use the scheme proposed by Entin-Wohlman et al [8| for U oo. 
The choice is justified by the simple analytical formulae which have an easy implementation. 

The triangular system is a good tool for studying the interplay between correlation and 
interference processes. The non-interacting conductance (curve no. 3 in Fig.l) shows a typical 
Fano line presenting both a Fano zero and a perfect constructive interference. This is the 
result of the interference between the partial waves passing the dot and the reference arm. 
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FIG. 1: Conductance through the triangular interferometer (see inset): (1) using exact formula 
(Eq.8), (2) with Ng's ansatz (Eq.l2), and (3) for non-interacting case. The parameters are Ei = 
E2 = Epermi = 0, F = 0.025, ti2 = tid = t2d = 0.02, T = 10"^'' (measured in units of half-band 
width). The horizontal line represents the asymptotic value for E^ going to ±00 . 



In the Kondo regime the dot transmits through the Kondo peak and the transmission 
phase is "frozen" at 7r/2 (as measured also experimentally in jisl. fioll). The result consists 
in a much slower variation of the interference conditions giving rise to a Fano line of reduced 
amplitude compared to the non-interacting case. This is known as the Fano-Kondo effect 
and is described by the curve no.l in Fig.l. 

Significant differences can be noticed between the conductance obtained by using the 
exact formula eq.(8) and the Ng's approximation eq.(12) 2^. The differences are rather 
large, indicating that the Ng's ansatz does not capture well the combined effect of correlation 
and interference. Fig.l shows that the destructive interference is overestimated. When the 
dot is empty {E^, > 0.1) all the three curves coincide, as expected. The main control 
parameter (also in experiments) is the position on the energy scale of the atomic energy 
Ed which can be changed by applying an external bias . Asymptotically, for \Ed\ 00, 
the interacting site "d" is decoupled from the other sites; consequently the transport is 
performed only through the background branch (i.e., that one connecting the sites "1" and 
"2"). The background conductance is also plotted in Fig.l and represents the control limit 
of our calculations. In the range of the gate potential where the correlations are important, 
the temperature dependence is expected to be of Kondo- type ,i.e. the conductance increases 




FIG. 2: a) The crossing point exhibited by the current curves calculated at different temperatures 
in the range [T = 10^^*^ — 10^^]. Curve (1) is for the lowest temperature. Other parameters are 
the same as in Fig.l. The inset shows the same behavior for the T-shape system, b) The same 
curves calculated with Ng's approximation, which misses the crossing point. 

with decreasing temperature. However, when multiple paths are possible (and interference 
plays an important role), an increased transparency of the quantum dot may give rise, on 
the contrary, to a reduced conductance of the whole system. The interference conditions 
can be changed by a magnetic field (as in [g^) or by variation of the gate potential as we are 
doing here. 

The isotherms of the conductance g = dJ/dfi plotted in Fig.2 show two different temper- 
ature regimes separated by a crossing point: on the right the conductance decreases with T, 
and behaves oppositely on the left side. Our calculation determines the crossing point with 
an accuracy within numerical errors, and suggests that there is a gate potential such that 

dg{E^„T)/dT = 0. (13) 
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FIG. 3: The phase of GJJ^ (in units vr) for the same parameters as in Fig. 2; the stars indicate the 
gate potentials where the Kondo temperature equals the plot temperature. The curve (1) is still 
in Kondo regime at Ed = —0.5. 

The same crossing point was found for a more simple model, namely the T-shape system. 
The T-shape geometry consists of two coupled quantum dots, only one of them (the non- 
interacting dot) being connected to leads. The advantage of this more simple model is that 
the conductance is simply expressed by the density of states at the site connected to leads 
(details are found in, e.g. j^) and the effect of the interference is more obvious. Basically, 
the T-shape and the triangle describe the same Fano-Kondo physics; the triangle is however 
a technical challenge and the first step versus more realistic models. 

The DoS at the Kondo dot always shows the specific Kondo peak at the Fermi energy. 
However, the DoS at the coupling site may show a Kondo peak or dip depending on the 
constructive or destructive interference conditions, respectively. The two cases give opposite 
temperature behavior. The interference conditions gradually change with the applied gate, 
giving rise to regions with opposite temperature dependence separated by a crossing point 
(see inset of Fig. 2). The existence of this point seems to be a fingerprint of the Fano-Kondo 
effect in systems with interaction and interference. The use of the exact formula eq.8 is 
essential, since the Ng's approximation eq.l2 misses the crossing point (as can be noticed 
in Fig.2b). This shows that the Ng's approximations fails not only quantitatively but also 
qualitatively. 

If we intend to identify the gate interval where Kondo correlations are important, one 
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has to calculate the Kondo temperature. Another visualization of the Kondo region is to 
plot the phase of the dot Green function 

(f) = atan{ImGdd/ReGdd), (14) 

that is known to "freeze" at the value 7r/2 in the Kondo regime. The phase is plotted in Fig. 3 
for several temperatures, and on each curve, we mark the point where the Kondo temperature 
equals the temperature of the curve. The Kondo temperature was computed with the 
Haldane [s], [2l| formula. One can notice from Fig. 3 that the crossing point ~ —0.05 (see 
Fig. 2) is indeed in the Kondo regime. For higher temperatures we have checked that the 
isotherms no longer pass through the crossing point as the transport is no more governed 
by the Kondo physics. 

In conclusion, we have developed an exact current formula easy to apply to multiple- 
connected meso-systems. The advantage of the formula is that it uses only retarded Green 
functions and can be applied to any system with left-right and time reversal symmetry 
with one-site interaction. In the absence of the symmetries we give an alternative approach 
based on the Ng's ansatz, which nevertheless is not sufficiently accurate for describing the 
correlations. It has been found that the interplay between the Kondo correlation and the in- 
terference effect gives rise to different temperature behaviors depending on the gate potential 
applied on the interacting dot, separated by a crossing point. 
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